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The effective electric field experienced by the unpaired electron in the ground state of PbF, which 
is a potential candidate in the search of electron electric dipole moment due to some special char¬ 
acteristics, is calculated using Z-vector method in the coupled cluster single- and double- excitation 
approximation with four component Dirac spinor. This is an important quantity to set the upper 
bound limit of the electron electric dipole moment. Further, we have calculated molecular dipole 
moment and parallel magnetic hyperfine structure constant (Ay) of 20 'Pb in PbF to test the accu¬ 
racy of the wave function obtained in the Z-vector method. The outcome of our calculations clearly 
suggests that the core electrons have significant contribution to the “atom in compound (AIC)” 
properties. 


I. INTRODUCTION 

The well established model of the interaction of ele¬ 
mentary particles, the Standard Model (SM), is incom¬ 
plete as it cannot explain some of the well known phe¬ 
nomena of fundamental physics. One such phenomenon 
is the dominance of matter over antimatter in our uni¬ 
verse, although the SM treats matter and antimatter ex¬ 
actly in the same way [lj. The violation of two funda¬ 
mental symmetries: inversion symmetry (P) and charge 
conjugation (C), is one of the several conditions that can 
explain the matter antimatter asymmetry [sj. The CP 
violation within the SM originating from the complex 
quark mixing Kobayashi-Maskawa matrix is too weak to 
explain such an asymmetry. Therefore, the search for an 
extra CP violation (flavour-diagonal CP violation), which 
is absent in the SM is needed to explore new physics be¬ 
yond the conventional SM ffjj, 0] . Another fundamental 
symmetry is the time-reversal (T) symmetry, which is 
also violated with the violation of CP symmetry though 
a direct observation of the violation of T symmetry is yet 
to observe [f|. 

The electric dipole moment (EDM) of any elementary 
particle is a consequence of violation of both T and P 
as dipole moment is odd under P and even under T but 
spin is even under P and odd under T [f|. According 
to SM, the electron EDM is too small (less than 10 -38 ) 
to observe experimentally Q. Therefore, a measurable 
non zero EDM of electron would be the proof of an extra 
CP violation and the first direct observation of T viola¬ 
tion [§j. However, the intensive search for the electron 
EDM over the period of half a century have not drawn 
any conclusion to the final value of electron EDM, which 
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would have out-turned in the upper bound limit of the 
electron EDM for different quantum systems. Till date, 
the best limit of electron EDM in an atomic system is 
achieved from the T1 atom experiment (\d e \ < 1.6 x 10~ 2 ' 
e cm) Q. However, the discovery of Sandars 0 reveals 
that the effective internal electric field experienced by an 
electron is profoundly enhanced in heavy-atonr contain¬ 
ing polar molecule which makes these polar diatomics 
very promising candidate in the search of the P- and P,T 
-violating experiments and creates a dimension to ex¬ 
plore new physics. The latest best upper limit of electron 
EDM is set from the ThO experiment (|d e | < 8.7 x 10 -29 
e cm) by ACME collaboration 0 - This limit is one 
order lower in magnitude than the previous best limit 
(|d e | < 10.5 x 1CU 28 e cm), which is obtained from the 
molecular YbF experiment HH . 

Another advantage of using diatomic molecules in 
the search of electron EDM is the D doublet structure 
of the ground and metastable state of such polar di¬ 
atomic molecules, which gives an additional enhance¬ 
ment and due to this reason different molecules having 
2 Ei ground state (YbF flU - flfil ]. HgF 0 , HgH, BaF 
|l8j[). 2 ni ground state (PbF 0-0), 3 Ai metastable 
state (ThO OJUll, ThF+ HfH+ 0, PtH+ 

[Hj], WC [29l.l30jj)i et cetera have been proposed. Among 
these molecules, PbF has some interesting characteris¬ 
tics, which make it a strong candidate in the search of 
electron EDM. Pb is neither a lanthanide or actinide f 

element nor a transition d element and PbF has a 2 ni 

2 

ground state, which means the unpaired electron is in the 
7r orbital while in most of the other molecules, the un¬ 
paired electron is in the a orbital in their ground state. 
PbF, being a 2 ni ground state molecule, the spin angu¬ 
lar momentum of the unpaired electron contributing to 
the magnetic moment almost cancels the orbital angular 
momentum contribution to the magnetic moment. This 
leads to a smaller g-factor in the 2 ni state of PbF. The 
smaller g-factor (31„ 32j makes it very insensitive to the 
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background magnetic field and this leads to reduction in 
some systematic errors in the experimental observation 
of electron EDM [33] • The other molecules having small 
g-factor like PbF are ThO, ThF+, H£H+, PtH+, WC, et 
cetera but that is in their metastable 3 Ai state. On the 
other hand, 2 IIi state of PbF is a ground state, which 
is easy to synthesize experimentally as compared to the 
metastable states. The energy shifts of the levels of op¬ 
posite parity in the ground rotational state of 207 PbF due 
to the f 1 doubling is canceled by the magnetic hyperfine 
interaction as a repercussion the gap between two oppo¬ 
site parity levels is very small (almost degenerate) Il9j . 
Therefore, the molecule can be polarized very easily with 
the application of a weak electric field and opposite sign 
of the Q doublet component leads to the cancellation of 
some systematic error. 

The effective electric field ( E e ff ) experienced by the 
electron in an atom or a molecule, which is equally known 
as P,T-odd interaction constant (Wd = E e ff/\fl\) is a 
non-measurable quantity. On the other hand, it is very 
important to set the upper bound limit in the search of 
electric dipole moment of electron. Therefore, one has 
to rely on a very accurate theoretical method to calcu¬ 
late E e f f , precisely. The calculation of E e f / of a heavy 
diatomic molecule is not a trivial task as it requires simul¬ 
taneous inclusion of both the effect of special relativity 
and electron correlation due to the intertwined nature of 
these two effects. The relativistic coupled cluster method 
using four-component wavefunction meets these require¬ 
ments to fulfill the purpose [34j . There are two alternate 
choice to evaluate one electron response properties in the 
normal coupled cluster framework: (i) the expectation 
value approach and (ii) the derivative approach. The en¬ 
ergy derivative approach is superior to the other one as 
the calculated property value is very close to the prop¬ 
erty value calculated using full configuration-interaction 
(FCI) method. The Z-vector method is a technique in 
the derivative framework, which simplifies the complex¬ 
ity associated with it. It calculates the energy derivate 
in a size extensive manner and capable of rendering ac¬ 
curate wave function both in the near nuclear region and 
apart from the nucleus and its performance has already 
been tested [35j |. 

In this article, we have chosen Z-vector method in the 
coupled-cluster single- and double- excitation approxima¬ 
tion (CCSD) for the calculation of effective electric field, 
E e f f , experienced by the unpaired electron in the ground 
state of PbF molecule. The parallel magnetic hyperfine 
structure (HFS) constant of 207 Pb in PbF molecule is also 
calculated to judge the accuracy in the calculated E e ff 
values, since both of these properties need an accurate 
wavefunction in the near nuclear region. Further, we have 
calculated molecular dipole moment of PbF molecule and 
both the calculated HFS constant and molecular dipole 
moment are compared with the experimental value and 
all these results are compared with the values calculated 
by means of other theoretical methods. 

The manuscript is organized as follows. A brief 


overview of the Z-vector method including concise de¬ 
tails of E e f f and parallel component of the magnetic HFS 
constant are described in Sec. El Computational details 
are given in Sec. IIIII We presented our calculated results 
and discuss about those in Sec. m before making final 
remark in Sec. m Atomic unit is used consistently unless 
stated. 


II. THEORY 


The theoretical estimation of E e f / can be obtained by 
evaluating the following matrix element 

E eff = w d \n\ = (* a \£?^\* n }, (l) 

3 £ 

where f 1 is the projection of total angular momentum 
along the molecular axis and To is the electronic wave- 
function corresponding to Q state, n is the total number 
of electrons and H<j is the interaction Hamiltonian of d e 
with internal electric field and is given by 

H d = -2*cd e 7°7 5 p 2 , (2) 


where 7 are the usual Dirac matrices and p is the mo¬ 
mentum operator. 

The parallel magnetic HFS constant in Dirac theory 
can be derived by taking the z projection (along molec¬ 
ular axis) of the expectation value of the corresponding 
Hamiltonian and is given by 


A,, 


—* u 


on x n 


I To), 


( 3 ) 


where / is the nuclear spin quantum number, a are the 
usual Dirac matrices and jlk is the magnetic moment of 
the nucleus k. 

The dynamic part of the electron correlation effect is 
included by the exponential structure of the coupled clus¬ 
ter wavefunction, which is given as 

I T C c) = e T |$ 0 >, (4) 

where, |<f>o) is the ground state single determinant wave- 
function and T is the excitation operator. The form of 
the T operator is given by 


N 

t = t 1 + t 2 + ... + t jv = ^t„, ( 5 ) 


where 

T m = X] t il:: a l a \...a j a i , (6) 

a<b... 


where i,j(a,b) are hole (particle) index and are the 
cluster amplitudes corresponding to the cluster operator 
T m . The equations of cluster amplitude can be obtained 
by pre-projecting excited determinant with respect to 
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|<f> 0 ) °f the above equation. In coupled cluster single and 
double model, T is Tb + T 2 . The equations for Tx and 
T 2 are given as 

(<S>“\(H N e T ) c \$ 0 ) = 0 , ($%\(H N e T ) c \$ 0 ) = 0 , (7) 

where |<h“) and |<I>“ b ) are singly and doubly excited de¬ 
terminant, respectively and H^r is the normal ordered 
Dirac-Coulomb Hamiltonian. The subscript c means only 
the connected terms exist in the contraction between Hjv 
and T, which ensures the size-extensivity. As, the nor¬ 
mal coupled cluster (NCC) is not a variational method, 
the energy is not optimized with respect to the determi- 
nantal coefficients in the expansion of the many electron 
correlated wavefunction and the molecular orbital (MO) 
coefficients for a fixed nuclear geometry. Thus the calcu¬ 
lation of energy derivative in NCC framework needs to 
include the derivative of energy with respect to determi- 
nantal coefficients and MO coefficients in addition to the 
derivative of these two coefficients with respect to the 
external perturbation field. Thus one needs to calculate 
these terms for each external field of perturbation. How¬ 
ever, in Z-vector method, this can be avoided with the 
introduction of a de-excitation operator, A, where the 
A amplitude equations are linear and perturbation inde¬ 
pendent. The second quantized form of the A operator 
is given by 

N 

A = Ai + A 2 + ... + Ajv = A ra , (8) 

n 

where 

A m = ^2 X ab'.:. a Wj ■ ■ ■ a b a ai ( 9 ) 

a<b... 

where i,j(a,b) are the hole(particle) indices and X^ b ' are 
the cluster amplitudes corresponding to the cluster oper¬ 
ator A m . The detailed description of A operator and A 
amplitude equation can be found in Ref. [361 ]. In CCSD 
model, A becomes, A = Ai + A 2 . The explicit equations 
for the amplitudes of Ai and A 2 operators are 

(<h 0 |[A(7J w e T ) c ] c |$“) + ($ 0 |(^e T ) c |$“) = 0, (10) 

(<S> 0 mH N e T ) c } c \<S>%) + (<t> 0 \(H N e T ) c \Z?) 

($“|A|$“ b ) + (<S> 0 \(H N e T ) c \<S>%) = 0. (11) 

Although the term ($ 0 |(17Are T ) c |4>“)(4>“|A|4>“ b ) in A 2 
equation (Aquation fill) results into one disconnected di¬ 
agram but the diagram is not of the kind of closed with 
disconnected part, it is linked (for details see Ref. 1 H). 
This ensures the extensivity. The equation for energy 
derivative can be written as 

A E' = ($ 0 |(Oive T ) c |$o) + ($o|[A(Oive T )c] c |$o>,(12) 
where On is the normal ordered property operator. 


III. COMPUTATIONAL DETAILS 

The DIRAC10 [37] program package is used to con¬ 
struct one electron spinors, two-body matrix elements 
and one-electron property integrals. Gaussian charge dis¬ 
tribution is considered to take care of the finite size of the 
nucleus where the nuclear parameters [38 ] are taken as 
default value of DIRAC 10. Restricted kinetic balance 
(RKB) [39] is used to construct small component basis 
functions from large component basis. In RKB, the basis 
functions are represented in scalar basis and unphysi¬ 
cal solutions are removed by diagonalizing free particle 
Hamiltonian. The positive and negative energy solutions 
are generated in 1:1 manner by this formalism. We have 
done five different calculations (A-E) by varying basis 
function and number of correlated electrons. For Pb, 
dyall.cv3z [|(| and for F, cc-pCVTZ (4lJ basis is used 
and two different calculations are done by using 55 and 
73 number of correlated electrons and these are denoted 
by A and B, respectively. We have done three more cal¬ 
culations by using 55, 73 and 91 (all electron) correlated 
electrons where dyall.cv4z (40| and cc-pCVQZ [4l| are 
used for Pb and F, respectively and these calculations 
are denoted as C, D and E, respectively. The cutoff used 
for A, B, C, D and E calculations are 3500 a.u., 1000 
a.u., 70 a.u., 70 a.u., and 70 a.u., respectively. We have 
used the experimental bond length (3.89 a.u.) [HJ for 
the calculation of properties of PbF in its ground state. 


TABLE I. Dipole moment, parallel magnetic HFS of 207 Pb 
and effective electric field of PbF 


Basis 

H (D) 

A|| (MHz 

) 

E eff ( GV/cm) 

Z-vector 

Z-vector Expt. 

1311 Z-vector Expt. 

120, 31] 

A 

3.71 

9865 


36.6 

B 

3.72 

9962 


37.5 

C 

3.82 3.5T0.3 9968 

10147 

37.2 

D 

3.82 

10043 


37.9 

E 

3.83 

10121 


38.1 


IV. RESULTS AND DISCUSSIONS 

The properties described by Eqs. [TJ [3] are also known 
as “atom in compound (AIC)” properties as these proper¬ 
ties strongly depend on the electronic configuration of the 
give n (heavy) atom instead of the bonds between atoms 
[43l ]. In Table [T] we present the molecular dipole moment 
(/r), parallel component of magnetic hyperfine structure 
constant (A||) of 20 'Pb and effective electric field experi¬ 
enced by the unpaired electron of PbF. From Table Q] it 
is clear that our dipole moment values are in good agree¬ 
ment with the experimental value [3ll ]. We got values in 








4 


the range from 3.71 D (basis A) to 3.83 D (basis E) due 
to different basis and number of correlated electrons but 
this range fits well within experimental limit (3.5±0.3 D). 

The calculated parallel component of magnetic HFS 
constant of 207 Pb in PbF shows an excellent agreement 
with the experiment mm-, specially for E basis where 
the absolute difference between theory and experiment 
is only 26 MHz. The relative error of the parallel mag¬ 
netic HFS constant in five different calculations (A-E) is 
shown in Fig. [Q The highest and lowest deviations of 
Z-vector value from experiment are for basis A (2.86%) 
and basis E (0.26%), respectively. This trend in the de¬ 
viations (expressed in 8%) is expected. When we go from 
triple zeta (TZ) basis to quadruple zeta (QZ) basis with 
same number of correlated electrons (from A to C and B 
to D where number of correlated electrons are 55 and 73, 
respectively) the 8% decreases as QZ improves the con¬ 
figuration space more by including one higher angular 
momentum basis function than TZ. On the other hand, 
for same basis, if we include more electrons in correla¬ 
tion calculation (for TZ, from A to B and for QZ, from 
C to E), the 8% decreases as the more number of corre¬ 
lated electrons includes more orthogonal space to Dirac- 
Hartree-Fock space and thus includes more correlation 
contribution to the property value. It is also interesting 
to see that in TZ basis, when we go from A to B, the 
addition of 18 electrons (i.e., 4s+3d+4p core electrons of 
Pb) improves the parallel magnetic HFS constant by 97 
MHz. In QZ basis, as we go from C to D and D to E, the 
addition of 18 electrons (i.e., 4s+3d-|-4p and ls-3p core 
electrons of Pb, respectively) improves the A|| value by 
75 MHz and 78 MHz, respectively. From this observation 
we can conclude that the core electrons have significant 
role in the correlation contribution of parallel magnetic 
HFS value. 

In Table [TJ we present our Z-vector results of E e ff of 
PbF in five different calculations where the value ranges 
from 36.6 GV/cm to 38.1 GV/cm. We believe that the 
value in E basis (38.1 GV/cm) is the most reliable value 
of E e f f of PbF system as its corresponding parallel mag¬ 
netic HFS value has the closest agreement with experi¬ 



FIG. 1. Relative deviations between Z-vector and Expt. val¬ 
ues of parallel magnetic HFS values. 


ment. In E basis, the Z-vector magnetic HFS value has a 
uncertainty of 0.26%. So, considering basis set and other 
higher order correlation and relativistic effects, we can 
conclude that the E e ff of PbF is 38.1 GV/cm with 4% 
uncertainty. 

We compared our Z-vector result with other theoret¬ 
ically obtained values. From Table HU it is clear that 
our all electron value in QZ basis for both dipole mo¬ 
ment and parallel magnetic HFS of 207 Pb in PbF has the 
best agreement with experiment among all the other the¬ 
oretical values. Baklanov et al did two calculations with 
13 correlated electrons in spin-orbit direct Cl (SODCI) 
methods - one without outer core (OC) (44j correlation 
correction and the other with OC correlation correction 
[45j . The SODCI with OC correction (45| calculation 
gives better value for A|| of 207 Pb but gives poorer value 
of molecular dipole moment. It is worth to remember 
that Cl is not size extensive and thus does not scale 
properly with number of electrons. So, Cl is not a re¬ 
liable method for the system with a reasonable number 
of electrons, especially with heavy atom containing sys¬ 
tems. Recently, Skripnikov et al [211 1 have done two two- 
component (2c) coupled cluster calculations - one with 
single and double approximation (CCSD) and the other 
with CCSD with perturbative triples (CCSD(T)) correc¬ 
tion. In their calculations, Skripnikov et al have included 
only 31 correlated electrons and removed 60 electrons (ls- 
4f inner-core electrons of Pb) by using “valence” semilo¬ 
cal version of the GRECP scheme [4a, [47J. In valence 
semilocal version of GRECP, the components were con¬ 
structed for nodal valence pseudospinors by interpolat¬ 
ing the potential in the neighbor of pseudospinor node 
to avoid the singularity in the potential. The problem 
of valence GRECP approximation is that it can lead 
to “non-negligible” errors for valence electronic states 
due to the improper reproduction of nuclear screening 
[4(1 . Although the molecular GRECP calculations are 
two-component ones, the proper four-component wave 
function near the nucleus is restored at the nonvaria- 
tional restoration stage that can lead to small errors. 
On the other hand our all-electron calculations are four- 
component at all the stages of calculations. Although 


TABLE II. Comparison of molecular dipole moment, mag¬ 
netic HFS constant and E e ff of PbF 


Method 

p A,, 

(Debye) 

(2°7pb) 

(MHz) 

E e ff 

(GV/cm) 

SODCI(13e) [44] 

4.26 

9727 

33 

SODCI(13e)+OC [45] 

5.00 

10262 

37 

2c-CCSD(31e) [21] 

3.97 

10265 

41 

2c-CCSD(T)(31e) [21] 

3.87 

9942 

40 

4c-Z-vector(QZ, all electron) 

3.83 

10121 

38.1 

Experiment [20, jjl] 

3.5 ± 0.3 

10147 
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the authors in Ref. [ 2 TJ claim that the “contemporary 
full-electron studies have not yet been able to unambigu¬ 
ously surpass our approach when it comes to AIC and 
spectroscopic properties of interest”, we believe that the 
explicit treatment of core electrons is necessary for some 
AIC properties where the polarization of the inner core 
electrons plays an important role, which is evident from 
our calculated parallel magnetic HFS constant value. 

V. CONCLUSIONS 

In conclusion, we have applied Z-vector method in the 
coupled cluster framework to calculate E e f f experienced 
by the electron in the ground state of PbF molecule. The 
calculated molecular dipole moment and A|| of 207 Pb are 
in excellent agreement with the experimental values. As 
the calculated HFS constant is in very good agreement 
with experiment, we can say that our calculated E e ff 


= 38.1 GV/cm is most reliable as both require accurate 
wave function near the nucleus and expectation value of 
their operator are similar in structure. The core elec¬ 
trons have significant contribution in the calculated val¬ 
ues, which is evident from our calculated results. There¬ 
fore, it desirable to treat all the electrons explicitly to 
have much more accurate and reliable result. 
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